Temperature dependent resistivity of spin-split subbands in GaAs 2D hole system 
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We calculate the temperature dependent resistivity in spin-split subbands induced by the inversion 
asymmetry of the confining potential in GaAs 2D hole systems. By considering both temperature 
dependent multisubband screening of impurity disorder and hole-hole scattering we find that the 
strength of the metallic behavior depends on the symmetry of the confining potential (i.e., spin- 
splitting) over a large range of hole density. At low density above the metal-insulator transition 
we find that effective disorder reduces the enhancement of the metallic behavior induced by spin- 
splitting. Our theory is in good qualitative agreement with existing experiments. 
PACS Number : 73.40. Qv, 71.30.+h, 73.50.Bk, 72.80.Ey 



The observation of an apparent "metallic" behavior in 
low density two dimensional (2D) electron or hole semi- 
conductor systems has attracted a great deal of recent at- 
tention because of the long-standing belief that all 2D 
(disordered) electronic systems (at least in the noninter- 
acting limit) are insulators (i.e., infinite zero temperature 
resistance in the thermodynamic limit). Although the ex- 
perimental situation is controversial (to say the least) be- 
cause the "metallicity" of the 2D system must invariably 
be inferred from an extrapolation (to zero temperature) 
of finite temperature (usually T > WOmK) data, there 
is broad agreement on a number of intriguing and man- 
ifestly interesting features of the observed 2D resistivity, 
p(T,n), as a function of temperature (T) and 2D free 
carrier density (n), that require theoretical understand- 
ing irrespective of whether there is a true T — 2D metal 
or not. 

One of the most striking experimental anomalies in this 
subject is the strong temperature dependence of the resis- 
tivity of this effective 2D conducting phase at low carrier 
densities and in the relatively modest temperature range 
of O.IK — AK where the measured resistivity may change 
by as much as a factor of three whereas the corresponding 
resistivity change in this temperature range (the so-called 
Bloch-Griinciscn behavior) for a real 3D metal (e.g., Cu, 
Al) is miniscule (less than 1%). In fact, very weak tem- 
perature dependence of resistivity is also observed in 2D 
electron systems at higher densities, but at lower densi- 
ties (and in relatively clean samples) p(T) shows strik- 
ingly strong temperature dependence, not found in any 
known metallic systems. This very strong dependence 
of p(T, n) on T at low n, specifically in 2D p-GaAs hole 
systems, is the subject matter of this paper. 

There have been two alternate concrete physical mech- 
anisms proposed in the literature to explain the strong 
temperature dependence of p(T, n) in the 2D "metallic" 
systems. One mechanism Q is the possible temperature 
dependence of effective impurity disorder in the system, 
arising, for example, from the temperature dependence 
of finite wave vector 2D screening which could be very 
strong at low carrier densities. This mechanism is 



quite universal in the sense that all 2D carrier systems 
should manifest some aspects of temperature dependent 
screening since random (unintentional) charged impuri- 
ties (and often, particularly in modulation doped 2D sys- 
tems, remote dopants as well) are invariably present in 
semiconductors. We have earlier shown through detailed 
transport calculations Q that the temperature depen- 
dent disorder scattering could in fact provide a reason- 
able semi-quantitative explanation for the observed tem- 
perature dependence of 2D resistivity behavior in n-Si 
MOSFETs, p-SiGe, and p-GaAs 2D systems. The sec- 
ond mechanism for the temperature dependence of re- 
sistivity, which applies essentially only to 2D p-GaAs 
hole systems in a quantitatively significant manner, has 
been discussed extensively by the Princeton group and 
others This manifestly non- universal mechanism 

involves a temperature dependent transport mechanism 
arising from intersubband hole scattering between two 
spin-split hole bands, which are present in the GaAs va- 
lence band. The magnitude of this effect should there- 
fore correlate with the valence band spin-splitting as has 
been demonstrated experimentally. This proposal has 
been somewhat controversial, and there are experimen- 
tal claims (and counter-claims) for nonexistence (and 
the existence) of this spin-splitting induced intersub- 
band scattering mechanism. It is obvious that both of 
these physical mechanisms must be simultaneously oper- 
ational in the p-GaAs system because scattering by ran- 
dom charged impurities and valence band spin-splitting 
are both unavoidably present in real 2D GaAs hole struc- 
tures. What is therefore required is a (technically highly 
demanding) theory which includes both mechanisms on 
an equal footing, and in this paper we provide such a the- 
ory which, as we describe below, qualitatively explains all 
of the observed features of the temperature dependent 
p(T, n) of 2D GaAs hole systems. Our work convinc- 
ingly shows that in order to understand the 2D GaAs 
hole data comprehensively a theory such as ours includ- 
ing both mechanisms equivalently is necessary although 
one mechanism or the other may very well be dominant 
in specific experimental p-GaAs systems depending on 
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the details of disorder and spin-splitting (which would 
vary from sample to sample), thus resolving the exist- 
ing controversy of whether spin-splitting is a relevant 
mechanism in 2D "metallic" behavior or not. In all our 
discussions in this paper, consistent with the accepted 
(but, not quite rigorous) terminology in this subject, the 
word 2D "metal" or "metallic" behavior will only mean a 
strong temperature dependence of p(T) at a fixed density 
(above the so-called metal-insulator transition "critical" 
density), with dp/dT being positive at low temperatures. 

The valence band holes in GaAs-AlGaAs heterostruc- 
tures have large spin-orbit interaction due to the inver- 
sion asymmetry of the zincblende crystal structure and 
the real space asymmetry induced by the confining po- 
tential, producing substantial zero-magnetic-field spin- 
splitting The lack of inversion symmetry in the p- 
GaAs system can lift the spin degeneracy and produce 
two groups of spin-split 2D subbands of holes with dif- 
ferent band dispersions and transport properties. Re- 
cently, Papadakis et al. || have experimentally studied 
the relationship between the 2D metallic behavior and 
the spin-splitting of the p-GaAS quantum well systems 
by tuning the symmetry of the well confinement potential 
while keeping the total hole density constant. They find 
that over a large range of densities (in a regime where 
the samples exhibit metallic behavior) the symmetry of 
the quantum well plays an important role in the sense 
that the temperature dependence of p(T) seems to corre- 
late with the (well asymmetry tuned) spin-splitting. At 
high hole densities, far above the critical density for the 
apparent metal- insulator transition (MIT), the metallic 
behavior is very pronounced as the spin-splitting becomes 
large, and the two spin-split subbands are occupied with 
unequal densities. This enhancement of the metallic be- 
havior is induced by the temperature dependent inter- 
subband hole-hole scattering because at constant total 
density the intersubband scattering becomes stronger as 
the spin-splitting increases. Thus, the temperature de- 
pendent hole-hole scattering may, by itself, give rise to 
the apparent metallic behavior in the spin-split situation, 
at least for large carrier densities. In low density samples 
[^],[|, however, the spin-splitting becomes much smaller, 
and the change of the resistivity with temperature aris- 
ing from intersubband hole-hole scattering is suppressed 
even in the presence of large asymmetry in the confining 
potential. Strong metallic behavior of p(T) at low densi- 
ties cannot therefore be explained solely by spin-splitting 
effects. This indicates that additional scattering mech- 
anisms other than hole-hole scattering are required to 
explain the metallic behavior at low densities. Hamilton 
et al. H claim, based on the interpretation of their own 
2D GaAs hole data, that in the low density systems the 
metallic behavior is determined not by the spin-splitting 
induced intersubband carrier-carrier scattering, but by 
the magnitude of the low temperature resistivity, that is, 
by the effective disorder in the system. 



In this paper we investigate the 2D hole metallic be- 
havior observed in p-type GaAs systems including both 
disorder and spin-splitting effects. This system forms 
two spin-split hole subbands (without any band offset) 
with very different effective masses. Our goal is to cal- 
culate the temperature dependent resistivity in the spin- 
split two subband system. We use the Boltzmann trans- 
port equation approach to calculate the low tempera- 
ture ohmic resistivity of the two subband system tak- 
ing into account the single particle relaxation times of 
each subband (as arising from charged impurity scatter- 
ing), Ti, and the hole- hole intersubband scattering be- 
tween the two groups of holes, 77^, which relaxes the 
relative momentum of the two carriers ||,[l0) and there- 
fore contributes to the net resistivity of the two sub- 
band system. We calculate the single particle relaxation 
times Ti from a multisubband transport theory taking 
into account the long range Coulombic scattering by ran- 
dom static charged impurity centers invariably present 
in semiconductor structures. Since screening plays a cru- 
cial role in the temperature dependence of resistivity j^j 
we consider both hole-impurity and hole-hole Coulomb 
interactions being screened on an equal footing by the 
2D hole gas in the random phase approximation (RPA) 
using a multisubband screening formalism. In the multi- 
subband system screening is enhanced since all subbands 
contribute. The screening effect can, in the RPA, be in- 
cluded in terms of a matrix dielectric formalism |ll].[l2| , 

e ii,i'j'{l) = fiij$i'f - v (<l)Tkj(q)Fij,i'j'(q)> where i,j are 
the band indices, v(q) = 2ne 2 /(Kq) with the background 
lattice dielectric constant of the sample k, Fij^ji (q) is the 
subband form factor determined by the quantum confine- 
ment potential, and 11^ (q) are the polarizabilities cor- 
responding to the two spin-split (i,j = 1,2) subbands. 
Since suppression of screening by disorder scattering may 
be important (particularly at low densities) we include 
collisional broadening in the screening function through 
the Dingle temperature approximation [Q. 

We assume that the spin-split subbands have parabolic 
energy bands with spin degeneracy g s = 1 (in each band) 
and isotropic Fermi surfaces with effective masses mi and 
m,2- In our model the hole population of each subband, 
for a fixed total density n = n\ + 112, is determined by 
the effective mass ratio, i.e., n\jn-i = mi/m2, where rrii, 
Hi are the effective mass and population density for each 
group of carriers. Since there are two groups of holes we 
need to consider intersubband hole-hole scattering, which 
contributes to resistivity in addition to hole scattering 
by charged impurities. In particular, hole-hole scattering 
plays an important role when the carriers in each band 
have different masses and densities. In a two-subband 
system the resistivity without intersubband hole-hole 
scattering is given by p(T) = 1/ [<?i{T) + er 2 (T)], where 
a.; = nie 2 {Ti)/rrii is the individual conductivity of the 
ith subband. These Ti can be calculated by using the 
usual Boltzmann transport equation, which becomes a 
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set of coupled transport equations for distribution func- 
tions associated with each subband [Q. 

The resistivity in the presence of hole-hole intersub- 
band scattering of the two-band system is given by pjl0[ 



P{T) 
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where Ui is the conductivity associated with the zth 
subband with the scattering times r^, and a = 
e 2 nThhfn/mim2 is the conductivity associated solely with 
hole-hole scattering with an average mass given by to = 
(mini + m 2 n 2 )/n and a hole-hole relaxation time (jhh) 
for the relative momentum of the two spin-split car- 
rier systems. We now note an important feature of our 
spin-split two subband system: Spin in each subband is 
conserved and therefore, Coulomb scattering, being spin 
independent, cannot cause real intersubband scattering 
(which would change the spin index of the scattered car- 
rier). Thus, only intra- subband transitions contribute to 
the intersubband hole-hole relaxation time, and we find 
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where g s — 1 is the subband degeneracy and 
J(p) = p/(27r) d6 sin 6W(q)[f(p)] 2 , where f(p) = 
(1 + p 2 + 2pcos8) 1 ^ 2 with p = (ni/n 2 ) 1 ^ 2 and q = 
(47r/ g s )y / nin2 sin(0)//(p), and the collisional proba- 
bility W(q) = 2ir \v(q) / s(q)\ 2 . In long wavelength 
Thomas-Fermi screening approximation (TFA) we have 
W T FA(q) = (27r) 3 /Lg s (TO 1 + to 2 )] 2 in our 2D system. In 
the absence of any real hole transitions from one sub- 
band to the other the RPA dielectric function is given by 
e(q) = l — v(q)\n„ii(q) + X\22{q)\- The hole-hole scattering 
time Thh is symmetric with respect to the interchange of 
indices 1 and 2. Note that r^ 1 cx T 2 while t~ 1 oc T for 
T/Tf <C 1 but the two asymptotic regimes need not co- 
incide (additionally, the Dingle temperature suppresses 
temperature dependence for T -C To), and therefore 
the net temperature dependence could be quite complex. 
When Thh S> Ti, which happens at very low tempera- 
tures, the total conductivity becomes the sum of the con- 
ductivities of each subband. When hole- hole scattering is 
very strong (or impurity scattering very weak) , Thh -C Tj , 
Eq. (|l|) become p(T) — (n\oi + n\a{) / '(n 2 ! <7icr 2 ) ■ Thus, 
in these two limits Thh makes negligible contribution to 
the system resistivity. The resistivity has a linear tem- 
perature dependence in these limits (provided T -C Tp, 
and Dingle temperature effects are negligible). When 
Thh ~ Ti the temperature dependence of the resistivity 
is strongly affected by hole-hole scattering and becomes 
(for T < Tp) p T = Po + aT + bT 2 , where p = p(T = 0) 
and a,b are positive density dependent constants. Our 
numerical calculation shows a > b when > Thh an d 
a < b when Tj < Thh- 




FIG. 1. Calculated hole-hole scattering time Thh as a func- 
tion of the density ratio r = n-^jni for fixed total densities 
n = 1.2 x 10 n cm~ 2 , n = 2.5 x 10 10 cm" 2 at T = O.IK. Solid 
(dashed) lines are the results from the screened Coulomb po- 
tential within the RPA (TFA). Inset shows the spin-split en- 
ergy subbands. 

In Fig. 1 we show the calculated hole-hole scattering 
time as a function of the density ratio r = n 2 /ni for 
fixed total densities, n = ni + n 2 = 2.5 x 10 10 cto~ 2 , 
1.2 x 10 n cTO~ 2 . Inset shows schematically the spin-split 
energy subbands with different effective masses. We use 
a fixed effective mass to 2 = 0.3to g throughout this paper 
and create unequal population of the subbands (spin- 
splitting) by tuning mi > to 2 . (Since Eq. (g) is sym- 
metric with respect to the interchange of indices 1 and 2 
we show the results for n^jni < 1). Fig. 1 shows that 
the scattering time of RPA is different from that of TFA 
when the density difference is small. Since most experi- 
mental situations |5|(| correspond to 0.5 < n%jn\ < 1 the 
use of the proper screening function (i.e. RPA) is crucial 
(and long wavelength delta function scattering of TFA is 
inadequate), especially at low densities where the spin- 
splitting is very small. Fig. 1 also shows that at a given 
density and temperature the hole-hole scattering time 
increases as the population difference of the two bands 
increases, and the inverse scattering time r^ 1 diverges 
logarithmically as r — > 1. (This follows analytically from 
Eq. (||) also.) The scattering time has strong density 
dependence at a fixed temperature: the scattering time 
decreases as the density decreases. 

The effect of changing the confining potential is shown 
in Figs. 2 and 3, where we show the calculated frac- 
tional change in the resistivity p(T)/p(Q) as a function 
of temperature. In Fig. 2 the results for a high density 
sample n = 1.2 x 10 11 cto -2 are given. The density ratio r 
decides the strength of the spin-splitting. (Smaller r in- 
dicates larger spin-splitting.) For m = rt 2 the fractional 
change in the resistivity shows a linear T dependence and 
is determined entirely by the mobilities of each band. In 
this case we have exactly the same results as the single 
band model with g s = 2 and n — n\ + n 2 . When we tune 
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FIG. 2. Calculated fractional change in the resistivity 
p(T)/p(0) as a function of for n — 1.2 x W 11 cm~ 2 . Various 
curves correspond to r = 1.0, 0.83, 0.67, 0.59 (from bottom 
to top). The temperature dependent resistivity for a given 
density is shown in insets. 

mi (keeping mi constant) the two subbands have un- 
equal densities, An = n\ — ri2 = nil — + r). As 
r decreases the relative magnitude of the resistivity in- 
creases with temperature. Thus, the strength of the 
"metallicity" is enhanced when the two spin-split sub- 
bands have more unequal populations (i.e., larger spin 
splitting). This shows that hole-hole carrier scattering 
enhances the metallic behavior, particularly at high den- 
sities. (We get similar results for other high density sam- 
ples.) 

In Fig. 3 we show the results for a low density system, 
n = 2.5 x 10 10 cm~ 2 . At this density experiment [[| shows 
that the lowest temperature resistivity, po, is an order 
of magnitude larger than that at n = 1.2 X 10 n cm~ 2 . 
In addition, p increases by about 100% with increas- 
ing back gate bias even though the total carrier density 
is constant. (Note that in the high density samples po 
increases by < 40% with the same increase of the gate 
bias.) We believe that in low density systems the effec- 
tive disorder becomes larger when the confining potential 
becomes asymmetric (due to the back gate bias) . The im- 
portant role of disorder in the 2D metallic behavior under 
a back gate bias has also recently been emphasized in Si 
inversion layer experiments |l~3j | . To incorporate the 
disorder induced collisional broadening corrections we in- 
troduce the Dingle temperature Tjj in the screening func- 
tion. The pure RPA (Tjj = 0) case completely neglects 
collisional broadening effects on screening. The effect of 
the Dingle temperature is to increase the zero temper- 
ature resistivity (by suppressing screening for T <C To) 
as the splitting becomes large. Fig. 3 shows that the 
temperature dependence of the resistivity is strongly en- 
hanced when we consider a fixed Dingle temperature and 
hole-hole scattering. Without hole-hole scattering (and 
with a fixed Dingle temperature) the spin-split system 
still has a stronger metallic behavior than the correspond- 
ing one subband spin degenerate system (r = 1), but 
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FIG. 3. Same as in Fig. 2 for a low density systems with 
n = 2.5 x 10 10 cm' 2 . Solid lines are for the equal density case 
(ni = 712). Thin (thick) lines are the results without (with) 
hole-hole scattering. In (a) we use a fixed Dingle temperature 
(Tb = O.ITf) for different r, and in (b) we use different Dingle 
temperatures; T D = 0.1, 0.15, 0.2T F for r = 1.0, 0.83, 0.67, 
respectively. The insets show p(T) for the various case in the 
main figures. 

the strength of the temperature dependence is reduced 
compared to the results with hole-hole scattering [Fig. 
3(a)]. When we vary the Dingle temperature (i.e., in- 
crease the level broadening) the metallic behavior with- 
out hole-hole scattering decreases as the spin-splitting 
increases (Fig. 3(b)), which is precisely what is found 
experimentally [p|. Thus, without introducing the hole- 
hole scattering we can obtain qualitative agreement with 
experimental results for lower density samples including 
only disorder and temperature dependent screening. 

It is unclear whether our model including only spin- 
splitting induced hole-hole scattering and screening in- 
duced temperature dependent disorder scattering can en- 
tirely quantitatively describe the metallic behavior ob- 
served in low density 2D hole systems. But it is clear 
that our zeroth order model should be an essential part of 
any theory in this subject. We have established beyond 
any reasonable doubt that both spin-splitting induced 
hole-hole scattering and temperature dependent screen- 
ing play important roles in the observed "metallicity" of 
2D p-GaAs systems. 
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